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Abstract 

We obtain Gaussian upper bounds for heat kernels of higher order differential operators with 
Dirichlet boundary conditions on bounded domains in R'^. The bounds exhibit explicitly the 
nature of the spatial decay of the heat kernel close to the boundary as well as the long-time expo- 
nential decay implied by the spectral gap. We make no smoothness assumptions on our operator 
coefficients which we assume only to be bounded and measurable. 
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Introduction 

Off-diagonal Gaussian upper bounds for higher order differential operators with bounded measurable 
coefficients were first obtained by Davies j4]. The operators considered in that exposition were of 
order 2m on (R-^) with 2m > N. The higher order operators, expressed as: 

Hf{x) := (0.1) 

\<y.\<.m 
\P\<,n 

were shown to have heat kernels with off-diagonal bounds demonstrated in the inequality: 



, / N , « / \X — V\ 2-"-! \ 

\k{t,x,y)\ <cir^exp -C2- + cgt (0.2) 

Subsequently Barbatis and Davies p] were able to obtain optimal values for the constants C2 and C3 
in terms of the ellipticity ratio and dimension. 

In this paper we address the question of upper bounds on heat kernels generated by uniformly elliptic 
differential operators with Dirichlet boundary conditions on bounded regions of R^. We make the 
same assumptions of the coefficients as in [3], namely that they are measurable and bounded, and 
consequently find it more convenient to carry out the analysis with greater focus on the corresponding 
quadratic forms. We do however assume that the quadratic forms, on a bounded region C R^, 
satisfy the ellipticity condition: 

c-i||(-A)^/||^<Q(/)<c||(-A)^/||^ (0.3) 
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for each / G (ft) and c strictly positive. 

Although (|0.2p holds for operators on bounded regions with Dirichlet boundary conditions, the bounds 
do not reflect either the spatial decay near the boundary or the long-time asymptotics. We extract 
the manner in which 

k {t, x,y) as x,y dft 
and show that the heat kernel has off-diagonal bounds demonstrated in the inequality: 



\k (t, ^, y) I < ci(l - t-"^ dix, diy, exp l^^ca '"^ Jj - - stj 

where < ^ < m ^ ^ * is the spectral gap. 

The techniques we employ are close to those employed in [i], but we give more emphasis to the 
analysis of spatial derivatives of the heat kernel. Moreover the exponential time decay is deduced by 
exploiting the spectral gap. 

Sharp off diagonal heat kernel bounds were also obtained by Barbatis [2] in terms of a non-euclidean 
metric based on the coefficients of the operator, replacing the term \x — y\ in (|0.2p by d{x,y). In 
the case of bounded regions, an off diagonal bound was obtained for highly non-convex regions for 
the uniformly elliptic operator by Owen [6], in which he used the geodesic distance but boundary 
behaviour was not the focus of that analysis. 

Throughout this paper we will assume that c and Ci represent strictly positive constants. 
Notation 

Given a G representing the multi-index (ai, a2, as, ...) where 

\a\ :— ai + a2 + a3 + -I- 

we define the corresponding operator: 



d°'^xid"^x2d"^x3....d'^» XN 
and the set Va such that 

Va — {r : Ti < ai for all i} 
Moreover given any multi-index r in Va we define the vector factorial as: 

(:) - (::) (::) (::) ■■■ (::) 

The directional derivative of order m of an appropriately smooth function along a vector v in is 
expressed as 

^-J- := . • • V (V (V/.t>) .v) .w 
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Quadratic Form 



It is helpful to give an indicative though a non-rigorous formulation of the family of higher order 
operators that we focus on in this the paper. The operator is defined more completely through it's 
quadratic form. Given a bounded domain fl in we express the operator of order 2m > A'' as : 

Hf{x) := Yl {--^)^"'^D"{a^,pix)DPf{x)) (0.4) 

|o:|<m 
|/3|<m 

where aa^/s are complex bounded measurable functions. 
The associated quadratic form Q 

Qif) ■■= J2 [ a^,0{x)D!'f{x)D<-f{x) (0.5) 

\<x\<m''^ 
|/3|<m 

defined with domain equal to the Sobolev space W^"^ {Vl) will be assumed to satisfy the ellipticity 
condition with a strictly positive constant c 

c-i||(-A)^/||^<Q(/)<c||(-A)^/||2 (0.6) 
We define the the spectral gap 

s = mf — 

fec^i^) ll/lli 

Since we have made the asstunption that ^ < 1, it will be informative to track the dependency on 
this constraint by defining the quantity < e < 1 — ^ and 

7:=m(l-e)-f 

For a given point a; in we define it's distance from the boundary dfl 

dx ■= d {x , dil.) = inf |a; — y| 
yedQ 

We define the function g (t) such that 



t> 1 



t — s 



1 Boundary Behaviour 

Having imposed Dirichlet boundary conditions on our operator we expect the heat kernel k {t, x, y) 
to vanish at the boundary. The precise nature of this decay can be deduced by application of the 
Sobolev embedding theorem : 
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We consider the norm on Cq (fi) to be defined as 

iifii iimfii ^ \D''f{y)~D''f{x)\ 
ll/llcj 2^ 11^ /Hoc + sup K 

where n and k are the integer and the fractional parts of 7 respectively. 

Lemma 1.1 There is a strictly positive constant c such that for all f G VF™'^ (fl) and any unit vector 
V in 

|^|<;^rfSQ(/)'^ 11/11^ (1.8) 

for all X d 
Proof 

By applying Fourier transform to we have for all a; in 



j,6R« I (y - a=) h 



and consequently 



hence for a positive 



< crff I (a. + IIP)"^ + 1/(0 M^C 

It then follows from Cauchy-Schwartz that 



Optimizing over /i to find that 



< 

e 



"i-^ll/lli 

completes the proof. 
□ 

We proceed to find an upper bound for the heat kernel by applying lemma [01 to ft '■— e^^*/ but 
first we need a more comprehensive upper bound for | ^ q^^^ \ by applying the Spectral Theorem. 
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Lemma 1.2 If ft is e ^*/ for some f in (fl) then 

Qift)<Ht) ml (1-9) 

Proof 

The inequality follows from 

Qift)<\\He-"'\\ ||e-«*|| ll/lli 

and an application of the Spectral Theorem. 
□ 

We can now combine lemmas 11.11 and 11.21 to yield our upper bound for the heat kernel. 

Lemma 1.3 The heat kernel k {t,x,y) generated by the differential operator H satisfies the inequali- 
ties 

\k{t,x,y)\ < c (1 - ^)"' dZ d; when t < § 

\kit,x,y)\ < c (1 - ^)"' e-^* d2 dl when t > f 

Proof 

If ft :— e^^*/ then from lemma [Ol we have 

Choosing v appropriately and integrating yields 

J k{t,x,u)f{u) d^u\ <f^d2 ffW^ 11/11^-^ ll/tlll 
o 

< t; dl~g{t)'^e-^'^ Wfh 
This inequality gives us a bound for the L2 norm of k {t, y, •) 

\\k{t,y,-)h<f.^dl~g{t)'^e-^'^ 
which can then be applied to obtain the required bounds 

\k{t,x,y)\ = \ j k{\,x,u)k{\,u,y) d^u\ 
n 

< ||fc(i,x,.)l|2||fc(|,y,-)l|2 

< ldZdyg{lY-'e-^'^ 

□ 
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2 Gaussian Bounds 



Gaussian bounds exhibiting boundary decay can be given by interpolation between the bounds found 
in lemma [T751 and those obtained in Davies [J.The drawback of this method however, is that the 
presence of the term cst does not imply the the long-time exponential decay that we expect and 
would like to show. 

Dirichlet boundary conditions imply a positive spectral gap and hence exponential time decay. We 
give bounds that more concisely exhibit this behaviour. 

One of the key features of hypothesis (|0.6p is that many required operator inequalities can be reduced 
to proving the corresponding inequalities for polynomial symbols by applying Fourier transforms. 

Barbatis and Davies pjj make a stronger assumption on the operator coefficients. They assume that 
for a strictly positive q 

X! "0,a/3Ca^/3 < ^ Oa/? (a^) < 9 ^ aO,a/3^aC/3 (2.10) 

|a|— m |ok|— m |q;|— m 

|/3|=m \p\=rn |/3|=m 

almost everywhere in and the non-negative coefficient matrix Aq = £^o,af3 satisfies 

((-A)"/, 5) = / ao^a^pD'^fW^d^'x (2.11) 
a 

for aU functions f,g e (n). 

They obtained the necessary estimates for the polyharmonic operator by way of Fourier transforms 

and an application of the polarization identity on the co-efficient matrix. We do not make this 
assumption here. 

2.1 Twisted Quadratic Form Inequality 

Given Xq in fl and a unit vector a in R^, we define the bounded function ipxg,a on fl 

ipxo,a.{x) {x - Xo,a) (2.12) 

Moreover given a real number A and dropping the subscripts on ipxo.a we define the multiplicative 
operator e'^'^ 

(e^^f) (x) :=e^'>(-)/(x) 

e^'^ is a homeomorphism on (ft). 

The twisted operator H\tf, is then defined as follows: 
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with 

Dom (iJAvO - [e'^'^f ■ f e Dom (H)} 
More importantly we define the twisted quadratic form, Qx^p as follows: 



a,/3 



= I (x) (e-^'^Z?" (e^^/)) (e^^i?'^ (e"^^/)) x 

°./3 o 



By Leibniz, for each / G {Q) we can expand 



^«gAV.j^ ^ (2.13) 



It is then possible to show that the difference per (A) :— (/) — Q (/) is 

where the summation runs over all r in Va and s in Vjg but where either |r| < |q;| or \s\ < |/3|. 
The terms where both r = ex and s = (3 are incorporated in Q (/) 

Since the coefficients of the operator are uniformly bounded we have for c > 
and applying Cauchy- Schwartz 

l<|Q;| = |/3|<m 

From lemma [5^41 we see that for e < 1 
ber(A)| < c J2 



l<p<m 
l<p<m 



It follows from lemma KTSl with p = 

ber(A)| < cie (1 + 0) Q (/) + cie'-^"' (l + ^se^"-!)'" A2'"||/||2 

and simplifying 

ber(A)| < ci£ (1 + 0) Q (/) + c,e'-'"^ ([1 + 0s] A)^'" ||/||2 
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Using 

Re [Q if) - Qxi, if)] =Qif)- Re [Qm, if)] < \Q if) - Qxf if) 

we have 

[1 - cie (1 + 9)] Q if) < ReQxi, if) + ciei-^™ ([1 + Os] A)'" (2.14) 

Moreover 

\ImQx4, if)\ < cie (1 + 0) Q (/) + c^e'-'"'- ([1 + ds] A)'" \\f\\l (2.15) 
We note from lemma IA.5I that these inequaUties are true for all positive 9 

2.2 Spectral Gap 

From the outset we know that H has a positive definite least eigenvalue s. We define 

H - H - s 

similarly 

Hx^f e-^^^He^^f (2.16) 

and crucially 

QxAf) ■■=QM'if)-4f\\l 

It is easy to see that the real an imaginary parts of the newly defined twisted form, Qx^ satisfy 

[1 - cie (1 + 9)] Q if) < ReQxi, if) + de^-^" ([1 + 9s] A)'" ||/||2 (2.17) 

and 

|/mQ,^ (/)| < cie (1 + 0) Q (/) + c,e'-^"^ ([1 + Bs] A)^" \\f\\l (2.18) 

2.3 Twisted Semigroup Inequalities 

Lemma 2.1 There is a positive constant k such that for all c > k and all < p < I 

Q'x,^ if) Qa^ (./) + c (1 + p) (1 + s)'"- A^" 11/11^ 

is a sectorial form with vertex at and a semi-angle smaller then ^. 
Moreover 

j^Q if) < ReQxi, if) + c (1 +p) (1 + s)'™ A2"||/||2 
l+p 

Proof 

With 9 = '^^p^ and £ = 2^ the second part of the RHS of dTTSl) 

c,ei-2™([l + 0,]A)2"||/||2 

is less than 

c?"22™-i (1 + s)'" A2"||/||2 
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hence 

p\ImQ'^,i. if) I 



This then impUes that 



and consequently 
□ 

We define the operator H'^ 

where c > such that 

\Arg{H'^^^fJ)\<^ 

We recall the following corollary from Kato[5]. 

Corollary 2.2 Let (3 he angle such that j < (3 < ^ and set 

_ 1 
^ tan (3 

Then e~^^'^^ is an analytic semigroup in the sector 

Sp^{z: \Arg (z) [ < /?} 

moreover 

||e-^i.^^l| < 1 

Proof 

see[p492:[5]] 

□ 

From Davies [3] Theorem 2.38 it is evident that there is a positive constant Cp such that 

\\H',,^e-"'^^^^\\<-^ (2.19) 
Re z 

for all z £ Sp. 

In the following estimates we let Hx^ be defined by (|2.16p . 



< T^Q if) + pcr22"-i (1 + s)'"' X'"^\\f\\l 

= [1- cie (1 + e)] Qf + pcf"22™-i (1 + s)'" A2™||/||2 

< iieg^v^ (/) + (!+ P)(cr22™"'(l + =s)""A2™) 

= ReQx^ if) + C2 (1+p) (1 + s)'" A2™||/||2 
= ReQ'x^^if) 

\I^Q'x,^f \ ^ 1 
ReQ'x^^f -P 

\ArgQ',^^ if) \<t'^n-' G) f 
H'x^^ 0(1 + 5)'" A'™ 
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Lemma 2.3 There is a positive constant c such that 

\\e-">-^'\\ < exp [c (1 + s)'" A^™^] (2.20) 

Proof 

Let f e L"^ (il) and define ft := e"^^**/, then solving 

= -2ReQx^{ft) 

< 2c(l + s)^'"A2"||/t||2 

proves the claim. 
□ 



Lemma 2.4 Whenever P > 0, there is a positive constant ci such for any < a < 1 we have 

|jif^^e-^^**||+/3(l + s)'"A2™||e-^-^'*|| < ^e^(i+")(i+^)""^'"* (2.21) 



Proof 



Applying the triangle inequality gives us 

||i?A^e-^^^*||e-^(i+^)""^""* < ||i/;_^e^^^**|| +c(l + s)2'"A2"||e-^^^*| 
and it follows from (|2.19p and lemma [ 



Then the LHS of (P?^ is bounded above by 



^^e(i+,)-A^"t + c (1 + /?) (1 + s)'" A2™e^(i+^)""^""* (2.22) 



For positive real x we can re-write xe^* as 

1 

a 

for some a > 0. We can then observe that 



(l+a)xt _ —axt 



xe < 



^{l+a)xt g— 1 



a t 
and see that similarly 



~ at 



Substitution into (|2.22p attains the claimed inequality. 
□ 
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Lemma 2.5 Let f G L'^ (fi) and define ft e ^^'t^ f then there are positive constants ci and C2 
such that 

Qift) < ^e^^(^+^)""^""*||/||^ (2.23) 
at 

Proof 

From lemma [2?l1 we have that for some positive constants ci and C2 

Qift) < c,ReQ^^{ft) + C2{l + sf"'X'^\\M\l 

< c4H^^fth\\fth + C2{l + sf"'X'"^Ml 
By using lemma we can see that for < a < 1 



at 

< £lg2c(l+a)(l+s)^'"A^'"t||^||2 

at 



□ 



Corollary 2.6 There are positive constants ci and ci such that 



Proof 

By Substitution from 
and 

it can be seen that 
hence for < a < 1 



at 



Qift) = Qift) -sWftWl 



ft = e"-^^**/ = e~"^'*'e-^'f 



Qift) = e^^*Q(e-«^^'/) - s\\ft 



Qie-"'^'f) < e-^^' [Qift) + sWftWl 

< e-2^*[fLe-(i+")(i+^)^'"^='"*||/||2 + ,||;,||- 



Applying the estimate from lemma [2751 completes the proof. 
□ 
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2.4 Heat Kernel bounds 

Theorem 2.7 The integral kernel k\^^ (t,x,y) of e^^^^'^ satisfies the inequality 

for some positive constants c and ci. 
Proof 

For f E (ri) we define ft :— e^^^^^f. From lemma [TTT] 

l^^ll^^^d^.Qift)"^ \\ft\\2 



and from corollarv l2.6l 



1^1 < ^d:e(^)[-(i+^)^'"^^'"-2^1*||/||i-||/*ll2 

Recalling that /t = ^^Hx^tj: _ ^-st^-Hx^t j ^i^qj^ with lemma [^751 we have the estimate 

||/*||^<exp(ec2(l + s)2™A2™t)e-"*||/||^ 
Substituting this estimate into (|2.25p 

t 2 

Integrating along the path to the boundary 

I / k^.^{t,x,u)f{u) < ^^d2ef^=^'+^^""'""-^i1/||2 

n 

Following a similar argument to that in the proof to lemma [T751 

\k\,4, {t,x,y) I < ||A;a,0 (|,a;, •) hWkx,^ (|>y>-) II2 
to yield the upper bound 

□ 

Theorem 2.8 The integral kernel k{t,x,y) of e~^^ satisfies the inequality 

]kit,x,y) I <c,{l- ^)-H-^ dldle^v [-C. - 
for some positive constants ci and C2 and where s is the least eigenvalue and < 7 < m 
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Proof 

We demonstrate the proof in two stages. Firstly optimising over ■0 and then optimising over A 
From lemma [B .41 we know that 

recalling the definition of ■0 from (|2.12p , for some unit vector a 

i) [y) - tp {x) ^ (y-x,a) 

If a is such that 

'tp{y)-ip{x) ^ -\y-x\ 
then substitution into the estimate (|2.24p yields 

\k {t, x,y)\<-^ rf7rf7e[c.(i+.)^'"A---.]t-A|.-y| (2.26) 

Optimising the exponent of the RHS over A we find 

^2m-i _ 1^ ~ y\ 

~ 2toc2 (1 + s)^"* t 

Substituting back into (I2.26P we find there is a positive constant C2 

\k(t,x,y) I < _2-d7d7exp -C2 \'^-y\^-' ^ _ A (2.27) 

1, ((1 + .)-.)-^ ) 

Recalling that e = 1 — have the required estimate. 

□ 



Appendices 

A Polynomial and Operator Symbol Inequalities 

When perturbing our quadratic form Q, the resulting twisted form Q\^^ generates cross terms of 
of the form X^D"^. We can use ellipticity and Fourier transform to estimate these terms with the 
polyharmonic (— A) ^ for m > p + \ol\ 

wyn'^fwi < J (cii^rH+c2iArH) |/(^)|2rfA'i 

< ci(||(-A)^/||2 + l)+c2(l + A2™||/||^) 

However a more detailed decomposition of these polynomials is required to attain tighter bounds on 
j||2_ gQ f]^g^^ can show exponential decay of the heat kernel in long time asymptotics. 
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Lemma A.l If a , h , p and q are all positive constants then 



for all e > where Cp q := ( ] ^ [ ~1~ ) ' strictly positive and strictly less than 1. 



.1+1 
Proof 

Follows from maximising — ex^'^'^ and substituting x — ^. O 

Lemma A. 2 Given p > q > there is a strictly positive constant Cq^p such that for all f £ (fJ) 

||(-A)Vl|2<c,,,||(-Af /lb 

Proof 

Follows from the Spectral theorem. □ 

Lemma A. 3 If f & (f2), A G R and p is a positive integer then whenever r is a multi-index for 
which \r\ < p 

\lXP-\-\D-fh < s\\ (-A)5 + e-^|Ani/||2 

for any e > 
Proof 

From the isometry of the Fourier transform we see that 



then using \{i^Y\^ < l^pl''! and lemmaEH it follows that 

\\XP-\-\D-f\\l < 1 [e||pf + £-^|A|^^]|/(Opd^i 



< e\\i-A)i f\\l + e'^\\\'^\\f\\l 



Rescaling e recovers the required inequality. 
□ 
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Lemma A. 4 Ifr and s are two multi-indicies such that \s\ <p—l and \r\ < p where p is a positive 
integer and A G R then 

\\XP-M (_A)^ IIA^'-I^I (-A)^ /II2 < £|1 (-A)^ fWl + 22f-i£i-2f X^P \\f\\l 

for all e < 2 

Proof 

When |r I = p 
We apply lemma 

IIAP-I'^I (-A)^ /||2||A''-I^I (-A)^ fh < II (-A)« /II2 (e|| (-A)5 fh + e^^lAni/lb) 

= 6|| (-A)« fWl + [e-^ \\f Wfh II (-A)5 /II2) 
Using the estimate |a6| < j-a^ + eb'^ we show that the above is 

< ell (-A)5 fWl + l.e-^\X\^'P\\f\\l + e\\ (-A)5 /||^ 

< 2e||(-A)5/||^ + e-f^|A|^P||/||^ 

_ £±l£l 

Imposing the condition e < 1 and maximizing e p^i" over s we have the inequality 

IIAP-IH (_A)^ /||2||A^'-I«I (-A)^ /II2 < 2e|| (-A)§ /||i + £^-'1ApP||/||i 
The proof for this case is completed on rescaling e. 

When \r \ < p 

Follows from lemma ETSl 
□ 

Having completed this decomposition, we obtain an estimate for lower order operators in terms of the 
higher order operator. 

Lemma A. 5 There is a positive constant Ci such for all p > 0, 9 > and positive integer p < m 
II (-A)5 fWl + pX'^fWl <ci{l + 9)Q if) + CP i^l + -jj X'"^\\f\\l 

for all f G (ri) and with s equal to the bottom eigenvalue. 
Proof 

We can show by an application of lemma lA.ll that for each p > Q this is 

II (-A)5 fWl + pX'^fWl < II (-A)« fWl + pMI + P/i-"^ A2'"||/||2 (A28) 
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lemisort Applying lemma IA.2I and cUipticity whilst recalling that we have positive-definite spectral 
gap s 

II (-A)5 fWl + pX'nml < cill (-A)^ fWl + pm^,\\f\\l + pm{l + f,r-' X'"^\\f\\l 

< c,{l + e)Q if) - c^esWfWl + pmpWfWl + pm (1 + A2™||/||2 

we can then set a = 

> mp 

□ 

B Davies' Twisted Operators and the Canonical Functional 
Calculus 

The two definitions are consistent 



Lemma B.l H and have the same spectrum. 
Proof 

Given a sequence of functions /„ where {H — z) /„ ^ we have {H\^ — z) {e~^^ fn) 
□ 



Lemma B.2 For each f (fl) and z in the resolvent set of H we have 

(z - H^^r' f ^ e''^ iz-H)-\^^f 



Proof 

Let B = (z — H\^)^^ then 



which is just 

and out statement follows. 
□ 



B (z - H^.^) f = f 
Be"^'^ (z-H) e^'t'f = f 



16 



Corollary B.3 For f & Co {R) we have a canonical functional calculus for the twisted operator given 
by 

f (Hx^) : = e-^^f (H) e^^ 

Lemma B.4 If f G Co {R) and kf {x,y) is the integral kernel of the operator f (H) then f {H\^) has 
integral kernel kf^x^ {x,y) where 

kf,xi, {x,y) = e-^^^-hf {x,y)e^^^y^ 

Proof 

Let V G L'^ (O) then 

j kf^x^{x,y) v(y)dy = {f{Hx^)v){x) 
n 

= {e-^^f{H)e^^v){x) 

= e-^'^^^^y fc/(a;,y)e^'^(s')w(y)dy 
n 

□ 
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